Reduced fidelity susceptibility in the one-dimensional transverse field Ising model 
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We study critical behaviors of the reduced fidelity susceptibility for two neighboring sites in 
the one-dimensional transverse field Ising model. It is found that the divergent behaviors of the 
susceptibility take the form of square of logarithm, in contrast with the global ground-state fidelity 
susceptibility which is power divergence. In order to perform a scaling analysis, we take the square 
root of the susceptibility and determine the scaling exponent analytically and the result is further 
confirmed by numerical calculations. 
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I. INTRODUCTION 

Within the last few years, some concepts and tools 
in the field of quantum information science [1], such as 
entanglement and fidelity, have been introduced in the 
study of quantum phase transitions (QPTs) [2], which are 
described traditionally in terms of order parameters and 
symmetry breaking. However, without a priori knowl- 
edge of the order parameter, when a system undergoes 
QPT, its ground state will change dramatically. Since the 
definition of fidelity is the overlap between two quantum 
states, it's natural to introduce it in the field of QPTs. 
The most studied fidelity is the global ground state fi- 
delity [3-21], which is an overlap between two ground 
states with a slight change of the parameters. Therefore, 
to avoid the arbitrariness of the small change in numerical 
computation, Zanardi et al. introduced the Riemannian 
metric tensor [8, 9] , while You et al. suggested the fidelity 
susceptibility [10], which is an effective tool to study the 
critical properties. 

So far, in most works, people have concerned with 
the global ground state fidelity. We note that as far 
as in Ref [22] the authors have introduced a reduced fi- 
delity which is defined as the trace of the product of two 
matrices. This definition is not rigorous in quantum- 
information theory, and they require one of the density 
matrices to be a pure state to agree with the rigorous 
mixed state fidelity [29]. While in this paper, the re- 
duced fidelity we used is defined rigorously in quantum- 
information theory and is studied in [14, 23-26]. It has 
been shown that the reduced fidelity susceptibility (RFS) 
is an efficient indicator of QPTs. 

In this paper, we study the scaling and critical be- 
haviors of the RFS of two neighboring sites in the fa- 
mous one-dimensional (ID) transverse field Ising model 
(TFIM) . We note that the quantum critical properties of 
some other quantities like concurrence [27] and von Neu- 
mann entropy [28], and the global ground state fidelity 



susceptibility [11] have been studied in this model. How- 
ever, the scaling and critical behaviors of them are quite 
different from the RFS for two neighboring sites . It's 
found that the scaling and critical behaviors of the RFS 
take the form of square of logarithm as (InN) 2 for finite- 
size system and (ln|A — 1|) 2 (N is the total number of 
spins and external parameter A — 1 is the critical point) 
in the thermodynamic limit, in contrast with the global 
fidelity susceptibility [11], which takes the form of N and 
1 — A| -1 , respectively. The different scaling and critical 
behaviors lead to a different finite-size scaling analysis, 
and for RFS we have to take its square root which is 
logarithmic divergent to carry out the scaling analysis. 

The paper is organized as follows. In Sec. II, we briefly 
review some concepts in fidelity and its susceptibility, and 
then we introduce a general formula of the RFS derived 
in [25]. In Sec. Ill, the scaling and critical behaviors are 
obtained analytically for both finite-size and thermody- 
namic limit, respectively. Then, using the square root of 
the RFS, we perform a finite-size scaling analysis. The 
scaling exponent is determined analytically, and is con- 
firmed numerically. 



II. REDUCED FIDELITY SUSCEPTIBILITY 

Firstly, let us review some concepts in fidelity suscepti- 
bility. The Hamiltonian of a quantum system undergoing 
QPTs can be written as 



H(h) = Hq + hHi, 



(1) 



where Hi is supposed to be the driving term with a con- 
trol parameter h. The fidelity between two RDMs of the 
ground state p = p (h) and p = p(h + 5) is defined as 
[29] 



F (h, S) = trvV^Sp 1 / 2 . 
The fidelity susceptibility is given by [4, 10] 
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therefore we have F ~ 1 — x<5 2 /2. 

In this paper, we study the RFS of two neighboring 
sites, and the RDM is block-diagonal as p — Q\ Q2 
due to the specific symmetries in the ID TFIM, and q\, 
Q2 are 2x2 matrices. In [25], we have derived general 
formulas of RFS for density matrix of block-diagonal as 
P = ®"=if?i> where Qi's are 2x2 matrices. However, 
in this work the density matrix satisfies det Qi ^ and 
tvQi 7^ (i = 1, 2), then we give the formula under such 
case as follows [25], 
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Xi 



4trft 



(tr^) 2 -4 det g'i 



[d h det (gj)Y 
det (ft) 



(4) 



paper we consider N is even, therefore M is half odd in- 
teger. As is well known, there is a critical point exactly 
at A c = 1 in the thermodynamic limit. The transition in 
this model is of order to disorder type due to the compe- 
tition between the Ising coupling and the external mag- 
netic field. For A < A c , the ground state of the system 
is paramagnetic, while for A > A c the strong Ising cou- 
pling introduces a magnetic long-range order of the order 
parameter (a x ) to the ground state. 

The RDM of two neighboring sites takes the form 
of [30] 



q' = dhQ, and \i 1S the 'susceptibility' for block 
Qi. From the above formula we can see that the RFS is 
only depends upon to the first derivatives of the matrix 
elements. 



Pi,i+1 — 



fu + z_ 

z_ u- 

w z + 

z + w, 



(6) 



III. SCALING AND CRITICAL BEHAVIORS OF 
RFS IN ID TFIM 



in the basis {| TT>, I U>, I U>, I IT)} with | |) and | |) 
spin up and down. The matrix elements are 



The Hamiltonian of the ID TFIM reads 



M 



H * = - E l x ^. 



3+1 + ° J- 



-M 



(5) 



where cr" (a = x, y, z) is a Pauli matrix at site j, A is the 
Ising coupling constant in units of the transverse field, 
and periodic boundary conditions (<t_m = °m) are as- 
sumed. The total spins are TV = 2M + 1, and in this 



u± = \ (1 ± 2(a z ) + {a&D) , 
z ± = \{{^af)±{ala\)). 



(7) 



where the translation invariance of the ID TFIM is con- 
sidered, and we set i = for convenience in what follows. 
From Eqs. (4) and (6), the RFS can be written explicitly 
with these matrix elements as 



J 



1 



4(u + +u-) 
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[d x u + - o\U—j + 4 (d\z-) H 



,o ,2 . (wd x w - z + d x z + Y 
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U+U- 



(8) 



The following discussions of critical behaviors of the RFS 
is based on this equation. 

The Hamiltonian could be diagonalized by using 
Jordan- Wigner, Fourier and Bogoliubov transformation 
sequently (see [2]), and the mean magnetization in the 
ground state is given by [31] 



M 



1 1-Acos^ 



q=-M 

where uj c f >q is the dispersion relation 



(9) 



ui<t, g = \Jl + A 2 - 2Acos(/) g , 
4> q = 2nq/N, 



(10) 



with q half odd integer. The neighboring two-point cor- 
relation functions are calculated as [32] 
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FIG. 1: RFS as a function of A for various system sizes. The 
positions of the peaks approach to the critical point at A = 1 
as increase with system sizes. 
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A cos (2</> g ) — cos cj> q 



(11) 



It's known that the most important themes in critical 
phenomena are scaling and universality. For a finite sys- 
tem, there are no singularities unless the ground state 
level crossing occurs. In the following, we will calculate 
the mean magnetization, the correlation functions and 
their first derivatives, and get the finite-size scaling and 
critical behaviors of RFS. 



A. Finite-size scaling behavior of RFS 

The numerical results of the RFS for a finite-size sys- 
tem are shown in Fig. 1, from which we can see that the 
peaks become sharper and sharper as the system size in- 
creases, and their locations approach to the critical point 
at the same time. It's expected that there will be singular 
at A c if N is infinite. Then to study the scaling behavior 
of RFS, we define its maximum over A as Xm with the 
location A m . For a finite-size system, A m ^ A c , and yet 
if N is large enough, A TO is so close to A c that it's a good 
approximation for us to use A m = A c in computing Xm- 
The summations in Eqs. (9) and (11) converge quickly 
as N increases, which means for large N, the summation 
can be replaced by the integral as 
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Then it's easy to evaluate the integrals and get 
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However, the first derivatives of the above quantities are 
divergent with respect to N at A c , and we have 
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The main contribution to the above expressions in the 
large- iV limit arises from the summation around the zero 
point of sin(</> q /2), since 



i_E — 

N^q\ 



1 



lniV, 



(15) 



where S represents the singular area. Immediately, with 
the help of Eqs. (4) and (11), we get 



d{ala\ 



dX 



A=l 



16 

3^2 



lniV. 



(16) 



The matrix elements and their derivatives could 
be evaluated by inserting the above results into 
Eq. (7). Here we should specifically consider d\Z- — 
d\ ((ffoof) — {oqOx)) /4, the divergent behaviors of 
d\(<jQaf} and dx^af) at A c are both (\nN)/n, while 
their difference is convergent and could be integrated as 
a constant: 



~dX 



A=l 



1 

3tt' 



(17) 



Then, with Eq. (8), we get the divergent form of Xm with 
respect to N as 



Xm (N)=*Ai (In TV + Cl ) 2 + c 2 , 



(18) 



where c±, Ci are constants that could not be determined 
analytically, and the coefficient 



-4, 



27tt 4 



144tH 



1024 



7T 2 (9tt 2 + 32) (3tt 2 - 32) + 4096 

0.1485. (19) 



To perform a finite-size scaling analysis, we should find 
the logarithmic divergent quantity, and the above result 
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FIG. 2: Square of the maximum RFS over A as a function of 
N. The red dots are numerical results and the slope of the 
solid line is \f~A~l. 
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FIG. 3: Finite-size scaling analysis is performed. The RFS as 
a function of system size and coupling, collapses on a single 
curve for various system sizes. 



suggests us to take the square root of Xm for large N to 
get the logarithmic divergent form as 



Xm (N) 



1/2 



1 A x InJV + const. 



(20) 



1 /2 

To confirm this, we compute \ m numerically, and the 
result is shown in Fig. 2, from which we can see that the 
slope is excellently consistent with \f~A\, while the linear- 
relation is obvious for not large TV because ci is small with 
numerical computation. However the maximum global 
ground state fidelity susceptibility S m is of a different 
scaling form as In S m = 2 In iV+const [11] , thus In S m but 

1/2 

not S m will exhibit similar scaling behavior with Xm at 
the critical point. 



B. Critical behavior of RFS in the thermodynamic 

limit 



In the thermodynamic limit the total spin N is infinite, 
then the summations can be replaced by the integrals 
(12) and are evaluated as: 



1-A 



K(k) 



7T 

>$<>?)=— -K(fc) 

1 



1 + A 

7T 

1 + A 



E(fc). 
E(fc). 



^°f)=3^[K(*) (A-1)(2A 2 + 1) 

-E(fc)(A + l)(2A 2 -l)], (21) 

where K(k) is the complete elliptic integral of the first 
kind and E(fe) is the elliptic integral of the second kind 
with k — 2\[Xj (1 + A). Therefore the critical behavior of 
the RFS is determined by the elliptical integrals. At A c , 
we have k c — 2\fX~ c j (1 + A c ) = 1, and E(l) = 1, while 



K(fc) is divergent as k — > 1, its asymptotic behavior is 

1 , 1 



K (k) ~ In 



In- 



II -A| 



(22) 



from which we have K(fc c ) (1 — A c ) = 0. Then, obviously, 
at the critical point, with the above analysis we get the 
same results with (13) in the previous subsection. 

However, their derivatives are singular at A = 1. We 
take (<7 Z ) for example, 



d(a z ) A+l 



E(fc) 



A 2 



1 



dX ttA v ' ttA(A + 1) 
where we have used the following relations 
dK(k) E(fc) K(fc) 



K(fc), (23) 



dk 
dE(k) 
dk 



k k ' 

E(k) K(fc) 



(1 - k 2 )k 



(24) 



Then as A approaches to A c , E(fc) converges quickly to 1, 
and with Eq. (22) we get 



d\ 



-im. 1 



7T |X — A] * 



(25) 



which is logarithmic divergent as A — * A c . The derivatives 
of the rest correlation functions are 



dX 
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dX 



im. 1 
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where the logarithmic divergent terms come from the el- 
liptic integral K(k). It's noticed that the coefficients of 
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the logarithmic terms are the same with those in the pre- 
vious subsection for a finite-size system. Then we get a 
similar divergent form of \ as 

x (\)~A 2 (ln-^-^+d^ + d 2 , 

where A 2 = A\, and d\, d 2 are constants that could not 
be determined analytically. 

As we have derived in the above, the divergent form 
of the RFS around the critical point is (In TV) 2 for finite 
system and (In 1 1 — A|) 2 in the thermodynamic limit. In 
order to make a finite size scaling analysis, it's feasible for 
us to use the square root of the RFS, which is logarithmic 
divergent, and according to the scaling ansatz [33], the 
square root of RFS is a function of N" (A — A r „), which 

behaves as x (A m , N) 1/2 - x (A, N) 1/2 ~ Q[N V (h -h m )], 
where Q (x) cx In (x) for large x. This function is uni- 
versal and does not depend on system size N. Since 
A 2 = A\, we determine the scaling exponent v = 1 which 
is the same with the concurrence [27] and the global 
ground state fidelity susceptibility [11], and is confirmed 
numerically shown in Fig. 3. 

IV. CONCLUSION 

In summary, we have studied the scaling and critical 
behaviors of the RFS for two neighboring sites in the 



ID TFIM. We found that the divergent behaviors of the 
RFS take the form of (IniV) 2 for finite-size system and 
(In j A — 1 1 ) in the thermodynamic limit, which are dis- 
tinct from other quantities like the global ground state 
fidelity susceptibility [11], concurrence [27] and entan- 
glement entropy [28]. Then we use the square root of 
the RFS which is logarithmic divergent to carry out a 
finite-size scaling analysis. The scaling exponent is de- 
termined analytically, and is confirmed numerically. It's 
shown that, the RFS undergoes singularity around the 
critical point, thus indicate that the RFS can be used to 
characterize the QPTs. 
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